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ON THE CONVERGENCE RATES OF GAUSS AND CLENSHAW-CURTIS 
QUADRATURE FOR FUNCTIONS OF LIMITED REGULARITY 

SHUHUANG XIANG AND FOLKMAR BORNEMANN 



Abstract. We study the optimal general rate of convergence of the n-point quad- 
rature rules of Gauss and Clenshaw-Curtis when applied to functions of limited 
regularity: if the Chebyshev coefficients decay at a rate 0{n^^^^) for some s > 0, 
^ ~~f Clenshaw-Curtis and Gauss quadrature inherit exactly this rate. The proof (for 

I ^ Gauss, if < s < 2, there is numerical evidence only) is based on work of Curtis, 

Johnson, Riess, and Rabinowitz from the early 1970s and on a refined estimate for 
^~^ Gauss quadrature applied to Chebyshev polynomials due to Petras (1995). The 

^^ convergence rate of both quadrature rules is up to one power of n better than 

^_^ polynomial best approximation; hence, the classical proof strategy that bounds the 

^^ error of a quadrature rule with positive weights by polynomial best approximation 

^~> is doomed to fail in establishing the optimal rate. 

-(— > 

2 1. Introduction 

I— I Though Clenshaw-Curtis and Gauss quadrature are classical topics in numerical 

analysis, it is quite hard to track down a theorem that would establish the optimal 
<^ rate of the error E„{f) of the n-point rules for functions / : [—1, 1] ^ ]R of limited 

l/~j regularity. Here, regularity is most conveniently measured^ by the exponent s > 

■^ of a decay rate a^ = 0{m~^~^) of the coefficients a„i of the expansion 

Cn m=0 

Z^ in terms of the Chebyshev polynomials T,„ (x) of the first kind of degree m; the 

,__l prime indicates that the first term is to be halved. We say that such a function / is 

^ of class X'^ and claim that the error of both quadrature rules inherits exactly this 

• ^H rate: 

X E„(/)=0(n— 1). (1) 

C^ As noted by Bornemann (2010, p. 893), the case s = 1 can be found explicitly in the 

classical literature (we denote by E^{f) the quadrature error of Clenshaw-Curtis 
and by E^(/) that of Gauss): if / e X\ 

• Riess and Johnson (1971/72) proved £^(/) = 0{n~^); 

• Davis and Rabinowitz (1984, §4.8) gave a sketch that E^{f) = 0{n^^). 

It is a fairly straightforward exercise, however, to extend the approach taken by 
these authors to the case of general s > 0: an approach that starts from the bound 

00 

\Enif)\^ E l«m| ■ |E„(T„,)|. (2) 



^Some ways to determine s are discussed in §2. 
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Figure i. Numerical evidence that n-point Gauss quadrature has an 0{n ' ^) 
error rate for integrating the functions fs{x) = |^ — 0.3|'' (left: s = 0.5, right: 
s = 1.5) on the interval (-1,1): £j?(/s) (dots), E|^(/s) (circles), c,«"""^ (solid line). 



By using aliasing of under-sampled trigonometric polynomials, Riess and Johnson 
(1971/72) and Curtis and Rabinowitz (1972) showed, for Clenshaw-Curtis and 
Gauss quadrature, that En{Tm) is, up to some remainder, periodic in m with a 
period of 0{n) and an average modulus of 0{n ). Hence, provided the remainder 
can effectively be controlled, one would read off the rate (1). If it were not for 
this proviso, the story could end here; but the precise state of affairs differs 
considerably: 

• For Clenshaw-Curtis quadrature, the remainder is a term of higher order, 
indeed; its effective control established by Riess and Johnson (1971/72) for 
s = 1 easily carries over to s > 0; see §3 of this paper. 

• For Gauss quadrature, the sketch given by Davis and Rabinowitz (1984, 
§4.8) neglects the remainder. Since it is not of strictly higher order, the 
remainder is much harder to control: aliasing holds asymptotically up to 
m = o{n^'^) only; for larger m, phase errors of order 0(1) enter. 

Accordingly, to rigorously deal with Gauss quadrature, we split (2) after the first 
0{n^'-^) terms; the tail is then easily estimated by the decay of the coefficients and 
a simple uniform bound of E„ (T,„); see §4. Using the estimate of the remainder 
given by Curtis and Rabinowitz (1972), we are able to prove the rate (1) up to a 
factor logn for s ^ 2, whereas the case < s < 2 yields a suboptimal 0{n^^'^'^) 
bound. Using a refinement of the Curtis-Rabinowitz estimate due to Petras (1995), 
Xiang (2012) has recently eliminated the logarithmic factor for s ^ 2 (there is, still, 
no improvement in the case < s < 2); see §5. 

Summarizing, we have proved (1) for all cases except for Gauss with < s < 2: 

Theorem. Iff e X^ , the error of n-point Clenshaw-Curtis quadrature and, for s ^ 2, also 
that of Gauss quadrature have the rate 0{n~^~^). For < s <2, the Gauss quadrature 
error is (at most) of size 0{n~^^^^). 

Numerical experiments with fg (x) = \x — 0.3|^, which is of class X'^ (see §2), and 
various < s < 2 (as in Fig. 1) has led us to the conjecture that Gauss quadrature 
enjoys the same 0{n~'^~^) error rate as Clenshaw-Curtis also for < s < 2 in 
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general. We remark that these experiments also show that the 0{n^'^ ) error rate 
cannot be improved for any of the two quadrature rules. 

Quadrature vs. best approximation. In his detailed study of the almost equal nu- 
merical performance of the quadrature rules of Gauss and Clenshaw-Curtis for 
functions of various regularity types, Trefethen (2008) proved a suboptimal 0(n^^) 
bound for functions / £ X*. In the Gauss case he based his rate estimate on the 
classical bound |£^(/)| ^ 4£2„+i(/) (see, e.g., Davis and Rabinowitz 1984, p. 333) 
where E* (/) denotes the error of best approximation by polynomials of degree n; 
if / G X^ this allows the straightforward estimate (see, e.g., Rivlin 1990, Thm. 3.3) 

CO 

K{f)^ E \am\=0{n-^). 

m=n+l 

In the case f{x) = \x\ (which is of class X^) the estimate is sharp, since it is known 
by a theorem of Bernstein that (Varga and Carpenter 1985, Eq. (1.18)) 

lim nE!(|x|) =0.2801694990... . 

Hence, Clenshaw-Curtis and Gauss quadrature converge with a rate that is typi- 
cally one pow^er of n better than the one of pol5momial best approximation. 

2. Functions of class X'* 

It is well known (see, e.g., Davis and Rabinowitz 1984, §4.8.1) that the Chebyshev 
coefficients am of f{x) are given by the Fourier coefficients of /(cos 6): 

2 /-i f{x)Tm{x) 2 f" 

«m = — / ^ - dr = — f (cos 6) COS m6 do. 

Asymptotic analysis of Fourier integrals can now be used to determine the decay 
rate of the a^: e.g., the function fs{x) = \x — ^\^ with — 1 < ^ < 1 and s > is of 
class X*^ since by the method of stationary phase (Olver 1974, §§3-ii-3-i3) 

fl„ = --Tm{^){l - 0"'^r(l + s) sin(7rs/2)m-'^-i +o{m-'-^) (m -^ 00). 

Alternatively (but often less sharp), decay estimates of Fourier coefficients based on 
the smoothness properties of / can be used; e.g., (Zygmund 1968, Thms. II.4.12): 

Let f be defined on [—1, 1]. If f {cos t) is k — 1-times differentiable with 
a piecewise k-th derivative of bounded variation, then f E X^. 

Since all derivatives of cos t exist and are bounded by the constant 1, the smooth- 
ness properties of /(cos t) can conveniently be inferred from those of /(x) (but not 
vice versa). In particular, if / itself is A: — 1-times differentiable with a piecewise 
A:-th derivative of bounded variation, we still get / G X . 

Remark. Denoting the total variation of that piecewise A:-th derivative of / by V, 
Trefethen (2012, Thm. 7.1) proved the explicit bound^ 

2V 

|«m| ^ rZT ("^ ^^ + 1); 



^We use Knuth's notation of the n-th falling factorial power: a" = a{a — 1) ■■■{« — n + 1). 
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using it, Xiang (2012) rendered the rate estimate (1) in the explicit form 

if n is sufficiently large (and, for Gauss quadrature, k ^ 2); an estimate that would 
asymptotically be, for f{x) = \x\, just a factor of 2 off the true state of affairs. 

3. Convergence rate of Clenshaw-Curtis quadrature 

Clenshaw-Curtis quadrature on [—1,1] is the interpolatory n-point quadrature 
rule that is derived from the nodes 

^/c = cos f ^— j-TTJ {k = \,...,n). (3) 

Now, it is well known that from Tm {x) = cos(m arccos x) one reads off aliasing due 
to undersampling, that is, with^ m = Ijin — 1) + 2r and — (n — 2) ^ 2r ^ w — 1 

Tmi^k) = T2|r|(xfc); 

w^hich implies, since Clenshaw-Curtis is exact for polynomials of degree n — 1, 

l^iTm) = l^{T2\r\)=l{T2\r\)- 

Here, I^if) denotes the quadrature formula as applied to / and !(/) the integral. 
Therefore, as m ^ n ^ 00, the quadrature error E^{T„i) satisfies 

With / e X'% that is, a„i = 0{m~^~^) for some s > 0, we follow the ideas of Riess 
and Johnson (1971/72, p. 347) in estimating 

00 

|En (/)! ^ E Kl ■ I£«(^'?)l = 0(Si) +0(S2), 



where 



„ ~ ^ l/|4r2-l| -1^ ^ r^i _.-2^ 

^^-g„^^„ (2;(n-l)+2.)-^^ - ^^ " ,?„ f^ ^ °^" ^ 



Because of 

00-1 00-1 

^ W^^\^^' I^ /5TT = as + 1), (4) 

r=-oo 1^' -^1 y=i I 

we immediately see that Si = 0(n~''~^); hence we obtain the rate estimate 

E^{f)=0{n-'-') (s>0), (5) 

which proves the theorem of §1 in the Clenshaw-Curtis case. 

^Note that we do not need, for both quadrature rules studied in this paper, to consider odd numbered 
Chebyshev polynomials: all their integrals and quadrature errors vanish because of symmetry. 
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4. Convergence rate of Gauss quadrature I 

As substitute for (3) there are asymptotic formulas for the nodes x/^ of n-point 
Gauss quadrature (the zeros of the Legendre pol5moinial of degree n): a classical 
one of Gatteschi (1956/1957) is, writing (pf- = (4A: — l)Tz/{in + 2) for short,^ 

Xk = cos((pk+lcot{(pk)n-^ + 0{k^^n''^)) (1 ^ A: ^ n/2). (6) 

Using this and an 0(n^^) bound on the weights, Curtis and Rabinowitz (1972, 
p. 211) proved that the error in integrating the Chebyshev polynomials is^ 

2 
(-1)^'-^ — - +0{m^/n^)+0{mlogn/n^) -n< r < n, 



E„ (Tm) 



7T 

(-1)^'— + 0{m^ / n^) + 0{m\ogn/ n^) r = ±n, 



iiln -^m = i{4:n + 2) + 2r with —n ^ r ^ n and ; ^ 0. This way, aliasing holds 
asymptotically for m = o{n^'-^) only; for larger m, phase errors of order 0(1) 
will render the estimate useless. Still, because of |Tm| ^ 1 on [—1,1] we get the 
uniform bound \E^{Tni)\ ^ 4. We now estimate E^(/) = E'„ + E" by splitting the 
Chebyshev expansion at an index of the order 0{n^^'^) with some < e < 1 to be 
chosen later Using the uniform bound of E^{Tm) we thus get the tail estimate 

CO / 00 -1 \ 

We are left with estimating the first 0(n^+'^) terms of the Chebyshev expansion: 

E'n=E l«2-?l ■ \E^iT2c,)\ = 0{S[)+0{S'2), 
q=n 

where 

~ ^ l/|4r^-l| , f ^ 



M \r\tn 0'(4« + 2) + 2r)^^+i j-'^ ^tin 0X4« + 2) + 2r)«+i 



1S+I' 



q=n q=n 



From (4) we immediately see that S[ = 0{n '^ ^). Likewise, we obtain 

[0(nP-«)^) 0<s<2, 
[O(logn) s^2. 



7«+is; 



4Curtis and Rabinowitz (1972, p. 208) stated this result with 0{n^^) instead of 0{k^'^n^^) — citing 
as source Abramowitz and Stegun (1965, p. 787), who had however misstated the result of Gatteschi 
(1956/1957): Gatteschi's term 0{k^^n^^) reduces to 0{n^^) only for those nodes Xj. that belong to a 
fixed interval in the interior of [—1,1]. However, the calculations of Curtis and Rabinowitz (1972) are 
fairly easy to fix: in the end, their estimate of Ejf ( T,„ ) turns out to be not affected at all. 

^Curtis and Rabinowitz (1972) stated the remainder in the form 0(1 /n) + 0(log n/n) f or m = 0{n); 
the explicit dependence on m given here follows from noting that the quantities hj of their paper scale 
with m/n: the first remainder term estimates a weighted sum of /;?, the second a weighted sum of \hj\. 



Xfc 
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Summarizing, the optimized choice e = 1/2 results in the rate estimate 

rQ(„-3s/2) 0<S<2, 

En (/) = { , (7) 

[0{n~'-^ log n) s;j2. 
which proves the theorem of §i in the Gauss case up to a factor logn. 

5. Convergence rate of Gauss quadrature II 

Xiang (2012) observed that we can get rid of the logarithmic factor rn (7) by 
using a refined estimate of Petras (1995, Thm. 1 and p. 199): upon replacing the 
bound in (8) by a later, sharper one also due to Gatteschi (1987), namely 

:fc=cos((^fc + icot(<^,c)(2M + l)-2 + 0(A:^3„-i)^ (1 ^ fc ^ n/2), (8) 

and by using some improved, individual estimates of the w^eights, Petras proved, 
within the range m = 0{n^), that 

^^ui'^^i" ^ +0(mVw'^)+0(m^log(w)/w4) \r\<n, 

\E^{Tm)\ = ; 

y +0(ot/m^) +0(OT*/n^) +0(OT^log(n)/n'') \r\ = n, 
where 2n ^m = y(4n + 2) + 2r with \r\ ^ n and ^ / = 0{n). Thus, we obtain 

„l+e 

E;= X: \a2^\-\E^iT2^)\=OiS[)+OiS[)+0{S'2), 

q=n 

where S\ = 0(m~^~^) is defined as in §4 and 

«: ;>/|4r^-l|/n jf j/n J_ 

' k |.fe (7(4« + 2) + 2ry+^ ^ f-[ ^tin (y(4« + 2) + 2r)«+i ^ n«+2 ' 

(j=n " q=n 



By 



,Eji?^-o('°e»). jEr-oK--). 



and, for < e < 1, 0{n^ ^ log") = o(l) we get S^ = 0{n ^ ^). Likewise 

{0{n^^-')'/n) 0<s<4, 

\o(logn/n) s^4. 

Summarizing, though the optimal choice e — 1/1 just reproduces (7) for < s < 2, 
it results, this time, in the rate estimate 

E^{f)=0{n-'-^) (s^2), (9) 

w^hich finally proves the Gauss case of the theorem of §1. 



Luigi Gatteschi (1923-2007) worked for nearly 60 years on the asymptotics of the zeros of special 
functions with a focus on explicit, useful error bounds; see Gautschi and Giordano (2008). 
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6. Open problems 

We leave the following open problems as challenges to the reader; their solution 
would require further, significant technical refinements of the methods used in 
this paper: to prove that, for / e X'', 

• the convergence rate is 0{n~^^^) for Gauss quadrature if < s < 2; 

• \E^{f)/E^{f)\ and its reciprocal stay uniformly bounded (cf. Fig. i). 

Acknowledgements. The authors thank Nick Trefethen for his continuing interest in 
this work and for his comments on some preliminary versions of the manuscript. 
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